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ABSTRACT 

Deep optical spectra of the high excitation planetary nebula NGC 1501 and its W04 cen- 
tral star are presented. A recombination line abundance analysis of the central star's emission- 
line spectrum yields He:C:0 mass fractions of 0.36:0.48:0.16, similar to those of PG1159 
stars. A detailed empirical analysis of the nebular collisionally excited line (CEL) and optical 
recombination line (ORL) spectrum is presented, together with fully three-dimensional pho- 
toionisation modelling of the nebula. We found very large ORL-CEL abundance discrepancy 
factors (ADFs) for 2+ (32) and Ne 2+ (33). The mean value of -5100 K for the T e derived 
from He I recombination lines ratios is 6000 K lower than the value of 1 1 100 K implied by 
the [O III] line ratio. This result indicates the existence of a second, low-temperature nebular 
component which could account for the observed ORL emission. Electron temperature fluc- 
tuations (t 2 ) cannot account for the high ADFs found from our optical spectra of this nebula. 

A three-dimensional photoionisation model of NGC 1501 was constructed using the 
photoionisation code MOCASSIN, based on our new spectroscopic data and using the three- 
dimensional electron density distribution determined from long-slit echellograms of the neb- 
ula by Ragazzoni et al. (2001). The central star ionising radiation field is approximated by a 
model atmosphere, calculated using the Tubingen NLTE Model Atmosphere Package (Rauch 
2003), for abundances typical of the W04 nucleus of NGC 1501 and PG1 159 stars. The neb- 
ular emission line spectrum was best reproduced using a central star model with effective 
temperature T c ff = HOkK and luminosity = 5000 L Q . The initial models showed higher 
degrees of ionisation of heavy elements than indicated by observations. We investigated the 
importance of the missing low-temperature dielectronic recombination rates for third-row el- 
ements and have estimated upper limits to their rate coefficients. 

Our single-phase, three-dimensional photoionisation model heavily under-predicts the 
optical recombination line emission. We conclude that the presence of a hydrogen-deficient, 
metal-rich component is necessary to explain the observed ORL spectrum of this object. The 
existence of such knots could also provide a softening of the radiation field, via the removal 
of ionising photons by absorption in the knots, thereby helping to alleviate the over-ionisation 
of the heavy elements in our models. 

Key words: ISM: abundances - planetary nebulae:individual: NGC 1501 - stars:Wolf-Rayet 
- atomic data 



1 INTRODUCTION 

The high excitation, northern planetary nebula (PN) NGC 1501 
has been the subject of several recent tomographic and spatio- 



kinematical studies (Sabbadin & Hamzaoglu 1982; Ragazzoni et al. 
2001; Sabbadin et al. 000a). The central star of NGC 1501 has 
been classified as OVI/WC4/W04 (Aller 1976; Tylenda et al. 1993; 
Gorny & Stasinska 1995; Crowther et al. 1998), and has also been 



2 Ercolano et al. 



the subject of several studies, particularly due to it being one of the 
few known pulsating PG1159-type PN nuclei (Bond & Ciardullo 
1993; Bond et al. 1996; Ciardullo & Bond 1996). A detailed study 
of the ionised nebular gas, has however, not been carried out to 
date, due to the lack of suitable spectroscopic observations. In this 
paper we present a spectroscopic analysis and three-dimensional 
photoionisation modelling of the ionised gas in NGC 1501. Based 
on new optical observations and using the three-dimensional pho- 
toionisation code MOCASSIN (Ercolano et al. 003a) and the den- 
sity distribution derived by Ragazzoni et al. (2001), we construct 
models in order to obtain better constraints on the nebular physical 
properties and elemental abundances. 



It is a well known and widely studied problem that optical 
recombination lines (ORLs) in gaseous nebulae yield higher em- 
pirical abundances than collisionally excited lines (CELs; for a re- 
cent review see Liu 2002); this is tightly correlated with Balmer 
jump (BJ) temperatures, T C (BJ), being systematically lower than 
those from the [O III] nebular to auroral line ratio (Liu & Danziger 
1993). One possible solution to the problem is that temperature and 
density fluctuations within the nebular gas lead to underestimated 
elemental abundances from standard CEL analyses (Peimbert 1967, 
1971). While this may be a valid explanation for objects with low 
abundance discrepancy factors (e.g. Peimbert et al. 2004), for cases 
where T C (BJ) is more than a factor of 2 lower than T c ([0 III]) the 
method is no longer applicable (Liu 2002). Moreover, the t 2 values 
derived from T C (BJ) and T c ([0 III]) for some less extreme nebulae 
are found to be too small to reconcile the ORL and CEL abundances 
(Liu 2002). The main problem encountered with the t 2 formulation 
is the fact that the T e -insensitive IR fine-structure lines also yield 
much lower ionic abundances than those implied by ORLs (e.g. Ru- 
bin 1968; Liu et al. 2000; Liu & et 001a; Liu et al. 001b; Tsamis 
et al. 2004). The presence of hydrogen-deficient metal-rich knots 
of high density embedded in the normal nebular gas has been pro- 
posed as an alternative to temperature fluctuations for the case of 
the PN NGC 6153 (Liu et al. 2000). Dual abundance photoionisa- 
tion models have been constructed for PNe NGC 6153 (Pequignot 
et al. 2003; Tylenda 2003), Ml-42 and M2-36 (Tylenda 2003) and 
for the hydrogen deficient knots of Abell 30 (Ercolano et al. 003b). 
As discussed later in this paper, temperature and density fluctua- 
tions in a chemically homogeneous medium are not sufficient to 
explain the ionic abundance discrepancy factors (ADF's) derived 
from our observations of NGC 1501. 



The new observational data is presented in Section 2, followed 
by a detailed empirical analysis of the ORL and CEL data, which is 
presented in Section 3. In that section we also carry out a study 
of the central star spectrum, providing identifications and mea- 
sured equivalent widths of the detected stellar emission lines, as 
well as an abundance analysis for the central star wind. The three- 
dimensional photoionisation modelling procedures are set out in 
Section 4, while Section 5 is concerned with the results of the mod- 
elling. After an investigation of the applicability of the t 2 formalism 
to our results, carried out in Section 6, there follows a discussion 
of how the presence of hydrogen-deficient, metal-rich knots could 
resolve the remaining discrepancies between ORL and CEL abun- 
dances as well as the over-ionisation of the nebular gas as predicted 
by the modelling. Our final conclusions are given in Section 7 of 
the paper. 



2 OBSERVATIONS 

In this section new observations of the optical spectrum of 
NGC 1501 are presented. The nebular analysis and the three- 
dimensional modelling which follow are all based on the new data. 
IUE spectra for the PN central star were also retrieved from the 
archive to provide some indication of the stellar energy distribu- 
tion. These are also briefly discussed in this section. 

2.1 Optical Observations of NGC 1501 

The observations were carried out in dark sky conditions on the 
night of 2003 August 1st, using the double-armed ISIS long slit 
spectrograph mounted on the 4.2 m William Herschel Telescope 
at the Observatorio del Roque de los Muchachos, on La Palma, 
Spain. The gratings used were the R316R in the red and the R600B 
in the blue, giving wavelength coverage from 3400-5100 A in the 
blue and 5100-8000 A in the red. The spectrograph slit was aligned 
North-South, and the slit width used was 0.75 arcsec (the slit length 
was 4 arcmin), giving a spectral resolution of 1.5 A FWHM in the 
blue and 2.8 A FWHM in the red. Five exposure were taken, of 
1800 s each. 

The data were reduced using standard LONG92 procedures 
in MIDAS. The 2D frames were bias-subtracted, flat-fielded, and 
cleaned of cosmic rays. They were wavelength-calibrated using 
exposures of a CuNe+CuAr calibration lamp. Spectra were flux- 
calibrated using wide slit (8") observations of the standard stars 
HZ44 and Feige 110. Unfortunately wide slit observations of 
NGC 1501 were not acquired and, as a result, absolute flux cali- 
bration of its central star could not be carried out. 

The top and bottom panels of Figure 1 show the extracted 
blue and red spectra ofNGC 1501 after integration over ~68 arcsec 
(338 pixels) along the slit. 

2.2 IUE Observations of NGC 1501 

Low- and high-dispersion ultra-violet spectra are available for 
NGC 1501 from the IUE archive. These have already been de- 
scribed by Feibelman (1998). Only the low-dispersion data is use- 
ful as the high-dispersion spectrum seems to be very underexposed. 
Moreover, as noticed by Feibelman (1998), the nebular contribution 
to the IUE data is very small, as indicated by the weakness of the 
C III] AA1906,1908 emission doublet, which is normally one of the 
major nebular emission features. Therefore, the IUE spectra could 
not be used for the nebular analysis carried out in this work. Never- 
theless, the retrieved IUE SWP28953 and LWP08948 spectra were 
useful in providing some constraint on the luminosity of the central 
star, as described in Section 4.4. 



3 EMPIRICAL ANALYSIS 
3.1 The Nebular Spectrum 

A full list of observed lines and their measured fluxes is presented 
in Table 1. All line fluxes, except those of the strongest lines, 
were measured using Gaussian line profile fitting. The fluxes of the 
strongest isolated lines were obtained by integration over the de- 
tected profile. Column 1 of Table 1 gives the observed wavelengths 
corrected for the heliocentric radial velocity of the nebula as de- 
termined from the Balmer lines. The observed fluxes are given in 
column 2. Column 3 lists the fluxes after correction for interstellar 
extinction, according to 7(A) = 10 c(H/3)/ (A) F(A), where /(A) is 
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Figure 1. The observed optical spectrum of the PN NGC 1501 integrated over the angular extent of the nebula (68 arcsec) and normalised such that 
F(H/3) = 100. 



the standard Galactic extinction law for a total-to-selective extinc- 
tion ratio of R = 3.1 (Howarth 1983) and c(H/3) is the logarithmic 
extinction at H/3 (see Section 3.3). The ionic identification, labo- 
ratory wavelength, multiplet number, the lower and upper spectral 
terms of the transition, and the statistical weights of the lower and 
upper levels, are given in columns 4-10, respectively. All fluxes are 
given relative to H/3, on a scale where R/3 — 100. 

3.2 The Central Star Spectrum 

The upper and lower panels of Figure 2 show the blue and red spec- 
tra measured at the central star position. In each panel the upper plot 
shows the unsubtracted spectrum and the lower shows the spectrum 
after subtraction of the nebular spectrum. The central star spectrum 
was integrated over 5.8 arsec (30 pixels) along the slit, and the neb- 
ular spectrum subtracted from it was taken an equal area 10 arcsec 
(50 pixels) away from the central position. The Figure shows that 
some emission lines have both a stellar and a nebular component, 
but most are only produced by one or the other. As expected, broad- 
ening of the stellar lines is evident, while the width of the detected 
nebular lines is much narrower. 

Our 2.5 hour exposure on the 14th magnitude central star of 
NGC 1501 with the 4.2-m WHT's ISIS spectrograph, led to a very 
high signal to noise spectrum. The stellar emission lines are rela- 
tively narrow for a star of its class, and so many lines could be iden- 
tified. Emission lines from four ionization stages of oxygen (O V - 
O VIII), together with numerous C IV and He II lines, were found to 



be present. We detected a very weak C III 5696 A emission feature, 
whose EW of 0.40 A was 1120 times weaker than that of the C IV 
5801,12 A doublet. Identifications and measured equivalent widths 
(EWs) for the stellar emission lines are listed in Table 2. Our EWs 
show quite good agreement with those measured for a subset of 
the lines by Stanghellini, Kaler & Shaw (1994), though we note 
that our ionic attributions for the lines at 3689, 3933, 4780, 4931, 
5291 and 5471 A differ from those of Stanghellini et al. Using the 
WO star classification criteria listed in Table 3 of Crowther et al. 
1998, the EWs listed in our Table 2 yield a spectral type of W04, 
in agreement with the subtype assigned by Crowther et al. to this 
central star. 

The airmass at the time of the observations was quite large 
(1.44-2.13), so that the effects of differential atmospheric disper- 
sion on the narrow (0.75 arcsec) slit needed for the nebular ORL 
analysis led to significant losses towards the blue end of the cen- 
tral star spectrum. For this reason, neither the spectral slope nor the 
absolute continuum fluxes of the central star of NGC 1501 can be 
derived from the optical observations presented in this paper. 

Despite the g-mode pulsations detected in the central star of 
NGC 1501 (Bond et al. 1996) the emission-line optical spectrum 
of this star (see Figure 2) points towards the presence of a thick 
wind - no stellar absorption lines were detected at all. Koesterke & 
Hamann (1997) derived a mass-loss-rate of 5.2-10~ 7 M^yi^ 1 and 
a terminal stellar wind velocity of 1800 kms -1 , while Feibelman 
(1998) estimated the terminal stellar wind velocity to be in the 3260 
to 3460 km s -1 range from its IUE low-resolution spectrum. Fol- 
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Table 1. Observed and dereddened relative line fluxes, on a scale where H/3 = 100. The integrated observed H(/3) flux from Cahn et al. (1992) was dereddened 
using c(H/3)= 1.0 to give an integrated dereddened flux of /(H/3) = 5.26 ■ 10 -11 ergs -1 cm~ 2 
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Table 1. -continued 
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3p' 4D 


2 


4 


4541.75 


1.04 


1.24 


Hell 


4541.59 


4.9 


4f+ 2F* 


9g+2G 


32 


* 


4563.02 


0.16 


0.19 


Mg I] 


4562.60 




3s2 IS 


3s3p 3P* 


1 


5 


4571.23 


0.19 


0.23 


Mg I] 


4571.10 




3s2 IS 


3s3p 3P* 


1 


3 


4610.34 


0.11: 


0.13 


on 


4610.20 


V92c 


3d 2D 


4fF2* 


4 


6 


4613.89 


0.11 


0.13 


Nil 


4613.87 


V5 


3s 3P* 


3p 3P 


3 


3 






* 


Oil 


4613.14 


V92b 


3d 2D 


4fF3* 


6 


6 




* 


* 


on 


4613.68 


V92b 


3d 2D 


4fF3* 


6 


8 


4621.17 


0.087 


0.10 


Nil 


4621.39 


V5 


3s 3P* 


3p 3P 


3 


1 


4634.66 


0.51 


0.58 


NIII 


4634.14 


V2 


3p2P* 


3d 2D 


2 


4 


4640.97 


1.18 


1.34 


NIII 


4640.64 


V2 


3p2P* 


3d 2D 


4 


6 


4647.17 


0.49 


0.56 


cm 


4647.42 


VI 


3s 3S 


3p 3P* 


3 


5 


4650.76 


0.80 


0.91 


cm 


4651.47 


VI 


3s 3S 


3p 3P* 


3 


1 


4657.12 


0.20 


0.22 


[Fe III] 


4658.10 


F3 


3d6 5D 


3d6 3F2 


9 


9 




* 


* 


CIV 


4658.30 




5 


6 






4661.37 


0.29 


0.45 


Oil 


4661.63 


VI 


3s 4P 


3p4D* 


4 


4 


4685.99 


36.50 


40.30 


Hell 


4685.68 


3.4 


3d+2D 


4f+ 2F* 


18 


32 


4697.10 


0.13 


0.14 


Oil 


4696.35 


VI 


3s 4P 


3p4D* 


6 


4 


4711.72 


2.44 


2.66 


[Ar IV] 


4711.37 


Fl 


3p3 4S* 


3p3 2D* 


4 


6 


4740.52 


1.87 


2.00 


[Ar IV] 


4740.17 


Fl 


3p3 4S* 


3p3 2D* 


4 


4 


4861.84 


100.0 


100.0 


H4 


4861.33 


H4 


2p+ 2P* 


4d+2D 


8 


32 


4881.24 


0.061 


0.060 


[Fe III] 


4881.11 


F2 


3d6 5D 


3d6 3H 


9 


9 


4922.45 


1.11 


1.07 


He I 


4921.93 


V48 


2p IP* 


4d ID 


3 


5 


4932.05 


0.37 


0.35 


[O III] 


4931.80 


Fl 


2p2 3P 


2p2 ID 


1 


5 


4959.47 


416.0 


393.6 


[O III] 


4958.91 


Fl 


2p2 3P 


2p2 ID 


3 


5 


5007.30 


1250. 


1151. 


[O III] 


5006.84 


Fl 


2p2 3P 


2p2 ID 


5 


5 


5199.15 


0.930 


0.77 


[NI] 


5199.84 


Fl 


2p3 4S* 


2p3 2D* 


4 


4 




* 


* 


[NI] 


5200.26 


Fl 


2p3 4S* 


2p3 2D* 


4 


6 


5412.04 


5.01 


3.68 


Hell 


5411.52 


4.7 


4f+ 2F* 


7g+2G 


32 


98 


5453.01 


0.21: 


0.15 


SII 


5453.83 


V6 


4s 4P 


4p 4D* 


6 


8 


5518.04 


1.09 


0.76 


[CI III] 


5517.66 


Fl 


2p3 4S* 


2p3 2D* 


4 


6 


5538.04 


1.04 


0.72 


[CI III] 


5537.60 


Fl 


2p3 4S* 


2p3 2D* 


4 


4 


5666.98 


0.19 


0.12 


Nil 


5666.63 


V3 


3s 3P* 


3p 3D 


3 


5 


5674.72 


0.23 


0.15 


Nil 


5676.02 


V3 


3s 3P* 


3p 3D 


1 


3 


5755.87 


0.39 


0.25 


[N II] 


5754.60 


F3 


2p2 ID 


2p2 IS 


5 


1 


5802.58 


0.71 


0.45 


CIV 


5801.51 


VI 


3s 2S 


3p 2P* 


2 


4 


5811.60 


0.72 


0.45 


CIV 


5812.14 


VI 


3s 2S 


3p 2P* 


2 


2 


5876.14 


18.70 


11.40 


He I 


5875.66 


Vll 


2p 3P* 


3d 3D 


9 


15 


6037.34 


0.53 


0.31 


Hell 


6036.70 


5.21 


5g+2G 


21h+ 2H* 


50 


* 


6074.27 


0.37 


0.21 


Hell 


6074.10 


5.20 


5g+2G 


20h+ 2H* 


50 


* 


6102.35 


0.76 


0.43 


[KIV] 


6101.83 


Fl 


3p4 3P 


3d4 ID 


5 


5 


6118.58 


0.52 


0.29 


Hell 


6118.20 


5.19 


5g+ 2G 


19h+ 2H* 


50 


* 


6235.04 


0.40 


0.21 


He II 


6233.80 


5.17 


5g+2G 


17h+ 2H* 


50 


* 


6300.45 


0.27 


0.14 


[OI] 


6300.34 


Fl 


2p4 3P 


2p4 ID 


5 


5 


6312.51 


2.85 


1.48 


[S III] 


6312.10 


F3 


2p2 ID 


2p2 IS 


5 


1 




* 


* 


He II 


6310.80 


5.16 


5g+2G 


16h+ 2H* 


50 


* 


6393.54 


0.13 


0.07 


[MnV] 


6393.60 




3d3 4F 


3d3 4P 


10 


6 


6407.79 


0.48 


0.24 


Hell 


6406.30 


5.15 


5g+ 2G 


15h+ 2H* 


50 


* 


6527.78 


0.60 


0.29 


Hell 


6527.11 


5.14 


5g+2G 


14h+ 2H* 


50 


* 


6548.69 


6.57 


3.16 


[Nil] 


6548.10 


Fl 


2p2 3P 


2p2 ID 


3 


5 


6563.26 


665.0 


318.3 


H3 


6562.77 


H3 


2p+ 2P* 


3d+2D 


8 


18 


6583.95 


20.90 


9.93 


[N II] 


6583.50 


Fl 


2p2 3P 


2p2 ID 


5 


5 
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Table 1. -continued 



-^obs 


F(A) 


/(A) 


Ion 


<^lab 


Mult 


Lower Term 


Upper Term 


9i 


92 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 



6678.73 


6.97 


3.21 


He I 


6678.16 


V46 


2p IP* 


3d ID 


3 


5 


6717.04 


3.90 


1.77 


rs III 


6716.44 


F2 


2p3 4S* 


2p3 2D* 


4 


6 


6731.42 


4.31 


1.95 


rs in 


6730.82 


F2 


2p3 4S* 


2p3 2D* 


4 


4 


6795.87 


0.28 


0.12 


TK IV1 


6795.00 


Fl 


3p4 3P 


3p4 ID 


3 


5 


6891.36 


0.84 


0.36 


He I 


6890.88 


5.12 


Sg+ ?G 

J 6 T - ' 


12h+ 2H 


50 


* 


7065 82 


6.12 


2.51 


He I 


7065 25 


V10 


2p 3P* 


3s 3S 


9 


3 


7136.34 


35.70 


14.34 


[Ar III] 


7135 80 


Fl 


3p4 3P 


3p4 ID 


5 


5 


7159.16 


0.11 


0.045 


He I 


7160.56 




3s 3S 


lOp 3P* 


3 


9 


7178.31 


1.11 


0.44 


Hell 


7177.50 


5.11 


5g+ 2G 


llh+ 2H* 


50 


* 


7231.51 


0.34 


0.13 


CII 


7231.32 


V3 


3p 2P* 


3d 2D 


2 


4 


7237.20 


0.92 


0.36 


CII 


7236.42 


V3 


3p 2P* 


3d 2D 


4 


6 




* 


* 


CII 


7237.17 


V3 


3p 2P* 


3d 2D 


4 


4 




* 


* 


[Ar IV] 


7237.26 


F2 


3p3 2D* 


3p3 2P* 


6 


4 


7262.88 


0.21 


0.08 


[Ar IV] 


7262.76 


F2 


3p3 2D* 


3p3 2P* 


4 


2 


7282.00 


1.42 


0.55 


He I 


7281.35 


V45 


2p IP* 


3s IS 


3 


1 


7320.36 


1.91 


0.73 


[Oil] 


7318.92 


F2 


2p3 2D* 


2p3 2P* 


6 


2 




* 


* 


[Oil] 


7319.99 


F2 


2p3 2D* 


2p3 2P* 


6 


4 


7330.69 


1.55 


0.59 


[Oil] 


7329.67 


F2 


2p3 2D* 


2p3 2P* 


4 


2 




* 


* 


[Oil] 


7330.73 


F2 


2p3 2D* 


2p3 2P* 


4 


4 




* 


* 


[Ar IV] 


7331.40 


F2 


3p3 2D* 


3p3 2P* 


6 


2 


7531.01 


0.72 


0.26 


[CI IV] 


7530.54 




3p2 3P 


3p2 ID 


3 


5 


7592.85 


1.28 


0.45 


Hell 


7592.74 


5.10 


5g+ 2G 


10h+ 2H 


5 


* 


7618.89 


0.78 


0.27 


N V 


7618.46 


V13 


7g+2G 


8h+2H 


* 


* 


7751.64 


9.84 


3.52 


[Ar III] 


7751.06 




3p4 3P 


3p4 ID 


3 


5 


8046.37 


1.70 


0.54 


[CI IV] 


8045.63 




3p2 3P 


3p2 ID 


5 


5 



lowing the method of Kingsburgh et al. (1995), recently reviewed 
by Drew et al. (2004), we estimate the terminal stellar wind ve- 
locity by measuring the half width at zero intensity (HWZI) of the 
C IV AA5801,12 doublet. Kingsburgh et al. (1995) showed that in 
the case of Sand 1, a SMC W03 star, the velocity corresponding to 
the HWZI of C IV AA5801,12 doublet, corrected for the 10.65 A 
doublet splitting, was consistent with the terminal wind velocity 
measured from the black absorption edge of the C IV AA1548,51 
resonance doublet, measured in a high resolution IU E spectrum. 
The validity of the method was later confirmed for several other ob- 
jects. We measured a full width at zero intensity (FWZI) of 103 A 
for the C IV AA5801,12 doublet, which, following the Kingsburgh 
et al. (1995) approach discussed above, yields a corrected HWZI of 
46.2 A and a terminal stellar wind velocity of 2390 km s _1 for the 
central star of NGC 1501. 



3.3 Reddening 

The logarithmic extinction at H/3, c(H/3) = log[/(H/3)/F(H/3)], 
was derived from the Balmer line ratios, Ha /H/3, H7 /H/3, H<5 /H/3 
and He/H/3. The measured ratios, together with the correspond- 
ing value of c(H/3) are given in Table 3. The theoretical line ratios 
used were those given by Storey & Hummer (1995) for T c = 1 1 kK 
and N c = 10 3 cm~ 3 and we adopted the Galactic reddening law of 
Howarth (1983). Given the uncertainty introduced by the different 
flux calibrations for the two arms of the spectrograph, we expect 
large error bars on our measured Ha/H/3 ratio and on the corre- 
sponding value of c(H/3). We finally adopted the value c (H/3) = 1.0 
as this best corrected the 2200 A interstellar absorption feature. 
This value, which is used to correct observed fluxes for interstel- 
lar extinction throughout this work, is within the range of values 



Table 3. Logarithmic extinction coefficient derived from the Balmer decre- 
ment 



Balmer decrement 


Ratio 


c(H/3) 




Theory 


Observed 




Ha/H/3 


2.837 


6.65 


1.15 


H7/H/3 


0.470 


0.359 


0.92 


H5/H/3 


0.260 


0.195 


0.69 


H^/H/3 


0.160 


0.106 


0.86 


Adopted 






1.0 



"Theoretical line ratios are from Storey & Hummer (1995). 



reported in the literature for this object - 0.96: Kaler (1976); 1.2: 
Cahn et al. (1992); 1.1: Stanghellini et al. (1994); 1.11: Ciardullo 
& Bond (1996); 1.05: Sabbadin et al. (000a). 

3.4 A recombination line abundance analysis for the W04 
central star 

The W04 central star of NGC 1501 has a pure emission line spec- 
trum and so we have used the recombination line method of Kings- 
burgh et al. (1995, KBS95) to derive the relative abundances of 
helium, carbon and oxygen ions in its wind. For the case of the 
LMC W03 star Sanduleak 2, Crowther et al. (2000) derived He, C 
and O abundances for its wind via spherical NLTE line-blanketed 
modelling, which showed good agreement with those derived from 
the recombination line analysis of KBS95, so we expect our re- 
combination line abundances for the nucleus of NGC 1501 to be 
reliable. 
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Table 2. Identifications and measured equivalent widths (EW) of emission features in the spectrum of the central star of NGC 1501. The EWs were derived 
after subtraction of the nebular spectrum. 



-^obs 

(A) 


EW 

(A) 


Ion 


-^lab 

(A) 


Lower Term 


Upper Term 


3566 


1.5 


C IV 


3566.8 


7 


15 


3620 


3.8 


C IV 


3607.3 


6p 


9d 




* 


O VI 


3615.6 


6p 


7s 


3689 


7.9 


C IV 


3689.7 


6 


9 


3723 


11.2 


C IV 


3720.1 


7 


14 




* 


O V 


3698-3725 


3p' 3 D 


3d' 3 D° 


3811 


24.5 


O VI 


3811.4 


3s 2 S 


3p 2 P° 


3834 


21.2 


O VI 


3834.4 


3s 2 S 


3p 2 P° 


3933 


11.0 


C IV 


3929.3 


7 


13 




* 


C IV 


3934.4 


5s 2 S 


6p 2 P° 




* 


O VI 


3936.6 


9 


13 


3976 


5.0 


? 








4026 


0.70 


He II 


4025.6 


4 


13 


4121 


6.2 


O V 


4119-4124 


3s' 3 P° 


3p' 3 D 


4153 


1.6 


O V 


4153.3 


3s' 3 P° 


3p' 3 D 


4161 


* 


C IV 


4160.2 


6p 


12d 


4199 


1.26 


He II 


4199.8 


4 


11 


4221 


5.0 


C IV 


4219.0 


6s 2 S 


8p 2 P° 


4229 


* 


C IV 


4229.1 


7 


12 


4265 


0.94 


? 








4298 


0.17 


7 








4340 


4.2 


He II 


4338.7 


4 


10 


4412 


0.23 


? 








4442 


6.8 


C IV 


4441.1 


5p 


6d 


4479 


0.14 


O V 


4479.5 


6g 3 G 


7h 3 H 


4499 


2.8 


O VI 


4499.0 


8 


10 


4520 


1.5 


O VI 


4521.2 


7 


9 


4541 


3.0 


He II 


4541.6 


4 


9 


4554 


4.0 


C IV 


4554.3 


6p 


8d 




* 


O V 


4554.5 


3p' X P 


3d' : D 


4607 


0.65 


C IV 


4604.6 


7p 


lid 


4659 


85 


C IV 


4658.7 


5 


6 


4686 


107 


He II 


4685.7 


3 


4 




* 


C IV 


4685.4 


6 


8 




* 


C IV 


4688.9 


7 


11 


4736 


0.60 


C IV 


4736.0 


6d 


8p 


4780 


9.8 


O VI 


4773.3 


7s 


8p 




* 


C IV 


4785.9 


5d 


6f 




* 


C IV 


4789.3 


6p 


8s 


4861 


9.1 


He II 


4859.3 


4 


8 




* 


C IV 


4860 


8 


16 


4931 


9.4 


O V 


4930.3 


6h 


7i 


4981 


0.66 


O V 


4977-4984 


6f 3 F° 


7g 3 G 


5018 


3.1 


C IV 


5017.8 


5p 2 P° 


6s 2 S 


5092 


13.0 


C IV 


5092.7 


8 


15 




* 


O VI 


5083.9 


7d 


8d 


5114 


* 


O V 


5114.1 


3s J S 


3p 1 P° 


5168 


0.81 


O VI 


5166.7 


10 


14 


5291 


10.4 


O VI 


5290.6 


7 


8 


5325-5358 


1.9 


C IV 


5335.2 


8p 


14d 




* 


C IV 


5356.9 


vp 


lOd 


5413 


15.1 


Hell 


5411.5 


4 


7 




* 


C IV 


5411.0 


8 


14 


5471 


15.8 


C IV 


5470.7 


7 


10 


5518-5538 


0.95 


C IV 


5518.8 


vp 


10s 




* 


C IV 


5530.5 


8s 


13p 


5591 


65.5 


O V 


5572-5604 


3p 3 P° 


3d 3 D 


5651 


0.58 










5667 


1.05 


O V 


5668.8 


9 


13 




* 


VII 


5669.3 


8 


9 
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Figure 2. Observed optical spectrum of the central star. Each panel shows the spectrum detected in the centre of the slit (upper line) and the same spectrum after 
subtraction of the nebular spectrum (lower line). Due to atmospheric dispersion effects on our narrow-slit spectrum, the slope of the stellar energy distribution 
is not reliable. 



To overcome the lack of a spectrophotometric calibration for 
our central star spectrum, we derived dereddened emission line 
fluxes by normalising the equivalent widths listed in Table 2 to a 
continuum flux slope corresponding to the model atmospheres dis- 
cussed in Sections 4.4 and 4.5. The T cfI =110 kK, log g=6.0 plane 
parallel NLTE Tubingen model discussed in Section 4.4 has a flux 
distribution between 4000 A and 7000 A which can be aproximated 
by Fa oc A~™, with n=3.93, while the spherically symmetric CM- 
Fgen model discussed in Section 4.5 has n=3.85. We adopted a 
spectral slope of n=3.90 and normalised it to Fa = 4.58xl0~ 14 
ergs cm' 2 s _1 A -1 at 5470 A, corresponding to dereddening by 
c(H/3) = 1.00 the continuum-level V magnitude of 14.44 measured 
by Stanghellini et al. (1994). Column 3 of Table 4 lists the resulting 
dereddened line fluxes used for the abundance analysis, while col- 
umn 4 lists the adopted line emission coefficients, Q (= hua rec ), 
for T e = 5 x 10 4 K and N e = 10 11 cm" 3 , taken from Table 14 of 
KBS95. Column 5 lists the ionic emission measures, I/Q, derived 
from each line, as well as mean I/Q values for each ion when more 
than one line was available. The C 4+ emission measures derived 
from the four C IV recombination lines used for the abundance 
analysis show excellent agreement with each other. Their mean I/Q 
value was used to predict the fluxes of the various C IV lines that 
blend with the He II lines used for our analysis. After correcting for 
these blends, the I/Q's derived from the He II lines showed excel- 
lent inter-agreement. 



Summing all ionic species, we find number ratios of 
C/He = 0.44, O/He = 0. 107 and (C+0)/He = 0.55 for the central star 
of NGC 1501. These ratios are not dissimilar to the number ratios 
of C/He = 0.48 and O/He = 0.054 derived by Barlow & Hummer 
(1982) from a recombination line analysis of the WOl PN central 
star Sanduleak 3. Expressed as mass fractions, we obtain He:C:0 = 
0.36:0.48:0.16 for the central star of NGC 1501, very similar to the 
values that have been derived from absorption-line analyses of pul- 
sating PG 1159-type objects, stars which have been postulated to 
be the descendants of PN WC/WO central stars (see Werner, Heber 
and Hunger, 1991). 



3.5 Nebular Empirical Abundance Determination 

A number of CELs, which are useful as nebular N e and T c diag- 
nostics, as well as for abundance determination, were detected in 
the spectra of NGC 1501. Using multi-level (>5) atomic models 
and solving the equations of statistical equilibrium, temperatures 
and densities were derived from these various CEL line ratios. The 
diagnostic line ratios and corresponding values of N e and T e are 
given in Table 5. 

The [S II] (A6716/A6731), [O II] (A3729/A3726), 
[CI III] (A5537/A5517) and [Ar IV] (A4740/A4711) diagnos- 
tic line ratios yield electron densities which are in good agreement 
with each other, giving an average density and standard error of 
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Table 2. -continued 



"^obs 


EW 


Ion 


^lab 


Lower Term 


Upper Term 


(A) 


(A) 




(A) 






5696 


0.40 


C III 


5695.9 


3p 1 P° 


3d*D 


5806 


448 


C IV 


5801.3,5812.0 


3s 2 S 


3p 2 P° 


5865 


7.0 


C IV 


5865.4 


8 


13 


5947 


1.46 


C IV 


5946.9 


9 


19 


6066 


0.88 


O VIII 


6068.2 


9 


10 




* 


viii 


6064.2 


11 


13 


6084 


0.34 


vii 


6085.1 


10 


12 


6120 


0.28 


He II 


6118.3 


5 


19 


6136 


0.36 


O V 


6131 


4s' 3 P° 


4p' 3 P 


6150-6170 


1.9 


C IV 


6149.9 


8s 


12p 




* 


Hell 


6170.7 


5 


18 


6200 


4.8 


O VI 


6200.8 


9 


11 




=i= 


O VI 


6198.3 


10 


13 


6235 


0.46 


He II 


6233.8 


5 


17 


6332 


0.85 


O V 


6330.1 


3p' : D 


3d' 1 F° 


6408 


8.4 


C IV 


6404.4 


7s 


9p 




=i= 


C IV 


6408.7 


9 


17 




* 


Hell 


6406.4 


5 


15 


6461 


4.3 


O V 


6460,66 


3p 3 D 


3d 3 F° 


6502 


1.35 


O V 


6500.2 


3p 3 D 


3d 3 F° 


6563 


51.4 


He II 


6560.1 


4 


6 




=i= 


C IV 


6559.5 


8 


12 


6624 


1.2 










6682 


2.4 


He II 


6683.2 


5 


13 


6747 


4.3 


O V 


6747.0 


9 


12 




=i= 


CIV 


6747.5 


9 


16 


7062 


53.0 


C IV 


7062.4 


7 


9 


7382 


4.3 


Civ 


7380.4 


6p 


7d 


7434 


2.2 


O V 


7422-7443 


4s' 3 S 


4p' 3 P° 


7723 


187 


C IV 


7726.4 


6 


7 




* 


C IV 


7735.9 


8 


11 




* 


O vi 


7717.1 


8 


9 


7862 


12.0 


C IV 


7860.7 


9 


14 


7924 


1.5 


O vii 


7926.4 


9 


10 


7949 


7.1 


Civ 


7947.2 


6d 


VP 



(1040±200) cm . A higher density (2400 crri) is derived from 
the [O II] (A3727/A7325) ratio. As a result of the higher critical 
density of the AA7320,30 doublet, the density derived from this 
diagnostic is weighted toward high density regions of the nebula. 

The [O III] nebular-to-auroral line ratio gives a temperature 
of 11.1 kK, which is slightly lower than the 12.6 kK given by 
the [N II] nebular-to-auroral ratio. Results from the photoionisa- 
tion modelling in Section 5 show that T c ([0 III]) is almost equal 
to T c ([0 II]) and T e ([N II]), and therefore the elevated tempera- 
tures measured from the [O II] and [N II] line ratios could be due 
to recombination excitation of the relevant CELs (Rubin 1986). 

We calculated the contribution of recombination excita- 
tion to the [O II] A7320,7330 and A3727,3729 lines us- 
ing Equation 2 from Liu et al (2000), and the relation 
I R (A3727,3729)=7.7 I R (A7320,7330) at 10 kK. Using the 2+ /H+ 
abundance derived from CELs (see Section 3.6), the corrected 
7325 A/3727 A line ratio is 0.0514, giving a revised temperature of 
13.7 kK, compared to the value of 16 200 K that is derived when no 
corrections are made. 

Using the formalism of Liu et al. (001b), we derived 
T C (BJ) = 9.4 kK from the ratio of the nebular continuum Balmer 
discontinuity at 3646 A to Hll A3770. We note that the signal- 



to-noise ratio in this region of the spectrum is rather poor (see 
Figure 1) and so the error on T C (BJ) is quite large. Also, scat- 
tered starlight affects the measured magnitude of the Balmer jump. 
However, by comparison of the flux of the scattered C IV A5800 
emission present in the nebular spectrum (Figure 1) with the equiv- 
alent width of the stellar C IV A5800 line, we estimate that scattered 
starlight accounts for only about 12% of the observed continuum, 
and so its effect is small. 

We also use the weak temperature dependence of some He I 
line ratios to derive T e . Table 5 shows that the He I line ratios 
yield much lower temperatures (4.1kK-5.1kK) than those derived 
from the CEL diagnostics (H.lkK-15.2kK). The observed relation 
T e ([0 III]) > T e (BJ) > T e (He I) is what would be expected in the 
presence of hydrogen-deficient inclusions within the nebula (Liu et 
al. 2001b). 



3. 5. 1 Collisionally Excited Lines ( CELs) 

The ionic abundances derived from optical CELs are listed in Ta- 
ble 6. They were derived using a constant temperature of 11.1 kK 
and a density of 1000 cm~ 3 . As discussed in Section 3.5, recom- 
bination excitation may contribute significantly to the flux in the 
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Table 4. Wind abundance analysis for the central star of NGC 1501. Column 1 lists the recombining ion, column 2 lists the transition and column 3 lists the 
dereddened line flux obtained by normalising the equivalent width listed in Table 2 to a stellar continuum that has a flux of F A =4.58x ~ 14 ergs cm~ 2 s _1 A -1 
at 5470 A and a spectral slope corresponding to F A oc A~ 3 - 90 (see text). Column 4 lists the line recombination coefficient Q for T e = 5xl0 4 K and 
N c = 10 11 cm~ 3 (see Kingsburgh, Barlow & Storey 1995), while column 5 lists the derived ionic emission measure, I/Q. Values appearing in brackets are 
where He II line fluxes have been corrected for the contribution made by blended C IV lines, where use has been made of the mean I/Q derived for C 4 + . 



Recombining Ion Transition I (ergs cm 2 s 1 ) Q (ergs cm 3 s x ) I/Q (cm 5 ) 



C 4 + 4658 (6-»5) 

5471 (KW7) 
5801,12 (3d 2 P° ^3p 2 S) 
7062 (9-t7) 
Mean C 4 + 

He 2 + 4686 (4->3) 

(C 4 + 4685 (8->6), 4689 (11— »7)) 
5412 (7->4) 

(C 4 + 5411 (14-»8)) 
6560 (6->4) 

(C 4 + 6560 (12-*8)) 
Mean He 2 + 

O s + 4930 (7->6) 

5591 (3d 3 D^3p 3 P°) 
Mean O s + 

O e + 5291 (8-+7) 

6200(11-^9, 13—^10) 
Mean 6 + 

7 + 6085 (12^10) 

O s + 6068 (10->9), 6064 (13-»11) 



Table 5. Plasma diagnostics 



Ion 


Lines 


Ratio 


Result 








logA r c (cm- 3 ) a 


[S II] 


7(A6731)//(A6716) 


1.10 


2.96 


[On] 


7(A3729)//(A3726) 


0.83 


3.00 


[O II] 


7(A3727)//(A7325) 


16.85 


3.38 


[CI III] 


J(A5537)/J(A5517) 


0.95 


3.22 


[Ar IV] 


7(A4740)//(A4711) 


0.75 


2.77 








T c (K) b 


[Nil] 


/(A6548 + A6584)//(A5754) 


52.60 


12600 


[O III] 


/(A4959 + A5007)//(A4363) 


148.7 


11100 


He I 


7(A5876)/J(A4471) 


3.03 


4100 


He I 


7(A6678)/7(A4471) 


0.86 


5100 


He I 


7(A6678)//(A7281) 


5.88 


4800 


He I 


7(A5876)//(A7281) 


20.84 


6500 


BJ/H11 




0.125 


9400 








±3000 



a Assuming T c = 11 100 K; 

b Assuming log Af e (cm -3 ) = 3.0. 



7.29xl0~ 12 6.65xl0~ 25 1.096xl0 13 

7.21X10" 13 6.06xl0~ 26 1.195X10 13 

1.626X10" 11 1.39xl0~ 24 1.169xl0 13 

8.96X10- 13 8.30X10- 26 1.080xl0 13 

1.135±0.05xl0 13 

8.96xl0~ 12 

-+6.43xl0" 12 2.38xl0~ 25 2.70xl0 13 

(2.53X10- 12 ) 2.23xl0~ 25 (1.135xl0 13 ) 

7.21X10- 12 

-^5.47xl0~ 12 2.20X10 -26 2.49xl0 13 

(1.74xl0~ 13 ) 1.53xl0~ 26 (1.135xl0 13 ) 

1.16X10- 12 

-^8.40xl0~ 13 3.28xl0~ 26 2.56xl0 13 

(3.20X10- 13 ) 2.82xl0~ 26 (1.135xl0 13 ) 

2.58±0.12xl0 13 

6.45X10- 13 8.53xl0~ 25 7.56xl0 n 

2.75xl0~ 12 4.33xl0~ 24 6.35X10 11 

6.96±0.60xl0 n 



5.42X10- 13 1.05X10" 24 5.16X10 11 

1.35xl0~ 13 3.06X10- 24 4.41X10 11 

4.77±0.37xl0 n 

1.028X10" 14 3.02X10" 25 3.40X10 10 

2.69X10- 14 1.72X10" 24 1.56xl0 1() 
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Table 6. Empirical ionic abundances of heavy elements derived from CELs. 



X i+ /H+ 


Lines 


N**+/ n h+ 


N+ 


[N ii] AA6548, 6584 


1.30(-6) 


o+ 


[O II] AA3726, 3729 


6.37(-6) 


o 2 + 


[O III] AA4959, 5007 


2.96(-4) 


Ne 2 + 


[Ne m]AA3868, 3967 


6.59(-5) 


S+ 


[S II] AA6716, 6731 


7.85(-8) 


S 2 + 


[S III] A6312 


2.22(-6) 


Cl 2 + 


[CI III] AA5517, 5537 


7.00(-8) 


Cl 3 + 


[CI IV] AA7531, 8046 


4.20(-8) 


Ar 2 + 


[Ar III] AA7135, 7751 


9.47(-7) 


Ar 3+ 


[Ariv] AA4711,4740 


3.09(-7) 


Fe 2 + 


[Fe III] A4658 


3.82(-8) 



[O II] A3727,3729 lines, and so the 0+/H+ abundance derived 
from them should be considered an upper limit. As the fraction of 
O in the form of + is very small (<2 per cent), the error this intro- 
duces into the total oxygen abundance determination is negligible. 

3.5.2 Optical Recombination Lines ( ORLs) 

The ionic abundances derived from the measured ORL spectrum of 
NGC 1501 are listed in Table 7. The ORL abundances were derived 
using the mean of the measured He I temperature of 5.1 kK, which 
is expected to be more representative of the temperature in the cold 
knots that give rise to most of the ORL flux rather than the [O III] 
temperature or the BJ temperature. The weak dependence of ORL 
abundances on physical conditions means that the difference be- 
tween abundances derived adopting the three temperatures are not 
large. 

The He + /H + abundance was derived from the intensities of 
the A4471, A5876 and A6678 He I lines, weighted 1:3:1 accord- 
ing to their approximate intensity ratios. According to the formulae 
given in Kingdon & Ferland (1995), collisional effects are negligi- 
ble at the temperature of 5.1 kK adopted here. The He 2+ /H + abun- 
dance was derived from the He II A4686 line. 

Heavy element ORLs are well detected in our spectrum. Our 
adopted ORL ionic abundances for C, N and O are the mean of 
the abundances derived from the selected individual lines for each 
ion. In Table 7 the ionic abundances that appear between square 
brackets were left out of the final averages for reasons that will 
be discussed next. The C II A5342 line is clearly too strong and 
must be contaminated by a blend. For this reason, this line is not 
listed in Table 7. The V3 C II A7231 line is a component of a 3p-3d 
transition, less reliable than the 3d-4f C II A4267 transition, and is 
therefore excluded from the average. For C 3+ we have preferred to 
use only the abundance derived from the V18 line at 4186.90 A, 
as this is a 4f-5g transition, hence more reliable than the VI line 
at 4651.47 A, which is a 3s-3p transition. Only one C IV line is 
detected in our spectrum and this is blended with [Fe III] A4658; 
we attributed 45% of the total emission to the C IV 5-6 transition 
(see discussion in Section 3.7), obtaining C 4+ /H + = 1.36 ■ 10~ 4 by 
number. For 2+ only the abundances derived from 3d-4f transi- 
tions were finally used for the average, with the exception of the 
V63a component at 4357.25 A which gives a high 2+ abundance 
and may be affected by blending (Tsamis et al. 2004, found that 
this line is 2.7 times too strong in the spectrum of NGC 3242). 



Table 7. Empirical ionic abundances for NGC 1501 from ORLs. The abun- 
dances listed between square brackets were not used in the calculation of 
the averages. See text for discussion. 



Ion 


A(A) 


Mult 


trans 


I 




He+ 


4471.50 


V14 




3.758 


7.29(-2) 




5876.66 


Vll 




11.398 


7.32(-2) 




6678.16 


V46 




3.215 


7.22(-2) 




Average 








7.29(-2) 


He 2 + 


4685.68 


3.4 




40.299 


3.03(-2) 


He/H 










1.03(-1) 


C 2 + 


4267.15 


V6 


3d-4f 


0.987 


8.25(-4) 




7231.32 


V3 


3p-3d 


0.131 


[3.02(-4)] 




Adopted 








8.25(-4) 


C 3+ 


4186.90 


V18 


4f-5g 


0.276 


2.56(-3) 




4651.47 


VI 


3s-3p 


1.462 


[7.00G-3)] 




Adopted 








2.56(-3) 


C 3+ 


4658.30 




5-6 


0.100 


1.36(-4) 


N 2 + 


403^ 08 
tUJU .uo 


V39a 




0.192 


1.83(-3) 




4621.39 


V5 




0.100 


3.29(-3) 




5666.63 


V3 




0.124 


1.02(-3) 




5676.02 


V3 




0.154 


2.85(-3) 




Average 








2.25(-3) 


o 2 + 


3882.29 


V11.V12 


3p-3d 


0.369 


[8.90(-3)] 




4285.69 


V78b 


3d-4f 


0.186 


7.75(-3) 




4307.23 


V53b 


3d-4f 


0.151 


1.14(-2) 




4357.25 


V63a 


3d-4f 


0.108 


[1.48(-2)] 




4477.90 


V88 


3d-4f 


0.1 14 


1.05(-2) 




4610.20 


V92c 


3d-4f 


0.133 


7.75(-3) 




4661.63 


VI 


3s-3p 


0.451 


[3.52(-3)] 




4696.35 


VI 


3s-3p 


0.138 


[1.16G-2)] 




Average 








9.35(-3) 


Ne 2 + 


3777.14 


VI 




0.278 


2.20(-3) 




4397.99 


V57b 




0.109 


[3.26(-3)] 




4409.30 


V55e 




0.113 


[1.76(-3)] 




4413.22 


V65 




0.157 


[6.80(-3)] 




4430.94 


V61a 




0.136 


[4.98(-3)] 




4457.05 


V61d 




0.124 


[1.30(-3)] 




Average 








2.20(-3) 


Mg 2 + 


4481.21 


V4 




0.064 


6.67(-5) 



In the case of Ne , as previously noted by Liu et al. (2000), the 
3d-4f coefficients are derived assuming the fine-structure levels of 
the ground term 2p 4 3 P2,i,o of Ne 2+ to be thermalised. However, 
given that the mean N c derived for NGC 1501 is more than an order 
of magnitude lower than the critical densities of the 3 Pi (2.0- 10 5 
cm ) and Po (2.9-10 cm ) levels of Ne , the assumption of 
Boltzmann equilibrium is probably unrealistic. At the lower densi- 
ties of NGC 1501 we expect the 3 Pi and 3 P levels to be under- 
populated and the 3 P2 level to be overpopulated. As the strongest 
3d-4f lines originate from the 3 P2 level, their recombination coef- 
ficient would be underestimated leading to correspondingly over- 
estimated abundances. From our observations the mean abundance 
derived from 3d-4f transitions is a factor 1.7 higher than the abun- 
dance yielded by the 3s-3p VI multiplet. In the light of the above 
discussion we adopt the Ne 2+ /H + abundance derived from Ne II 
multiplet VI only. 
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Table 8. ICFs and total elemental abundances by number with respect to H 
derived from CELs and ORLs. 



Ion ICF Abundance 



KB94 


MC2 


KB 94 


MC2 


from CELs 








N 60.4 


252 


7.85(-5) 


3.28(-4) 


O 1.27 


1.13 


3.85(-4) 


3.42(-4) 


Ne 1.30 


1.06 


8.57(-5) 


6.98(-5) 


S 2.73 


14.4 


6.29(-6) 


3.31(-5) 


CI 2.83 


1.13 


1.98(-7) 


1.27(-7) 


Ar 1.02 


1.03 


1.28(-6) 


1.29(-6) 


Fe 


155 




5.92(-6) 


from ORLs 








C 1.02 




3.52(-3) 




N 1.02 




2.29(-3) 




O 1.30 




1.22(-2) 




Ne 1.30 




2.86(-3) 





3.6 Summary of Nebular Empirical Abundances 

Total elemental abundances were first of all derived using the ion- 
isation correction factor scheme of Kingsburgh & Barlow (1994, 
KB94) (see columns 2 and 4 of Table 8). For oxygen, + /H + and 
2+ /H + abundances are available from CELs, and the unseen 3+ 
was corrected for using 

O _ / / Hc++He 2 + \ 2/3 / 0+ + Q 2 + \ 

R \ He + J ' { H+ ) ■ W 

The ionisation correction factor derived from this is 1.27. For 
oxygen ORLs, only 2+ /H + is available, and to derive a total abun- 
dance it was assumed that the ORL + /0 2+ ratio was the same as 
the value of 0.022 derived from CELs, giving a total ICF of 1.30. 
It is worth noting at this point that the ionic fractions in the ORL- 
and CEL-emitting regions may actually differ, but, while we cannot 
constrain the exact amount, we expect that only a small fraction of 
oxygen to be in the form of + , even in the coldest regions (e.g. 
Ercolano et al. 003b), and so the errors in the determination of O/H 
should be small. 

The nitrogen abundance from CELs is rather uncertain - only 
[N II] lines are seen, and N + accounts for only a small fraction of 
the total nitrogen abundance. KB94 give N/H = (N+/H+)-(0/0 + ), 
which gives an ionisation correction factor of 60.39. Uncertainties 
of the empirical determination of the ICF for N are discussed to- 
ward the end of this section. In the case of the ORLs, N 2+ and N 3+ 
abundances are available. The unseen N + is corrected for by as- 
suming that N/N + =0/0 + , which results in an ionisation correction 
factor of 1.02. 

No carbon CELs are seen as only optical nebular spectra were 
available. As discussed in Section 2.2, the IUE observations of this 
object are unsuitable for nebular analysis. However, nebular C II, 
C III and C IV ORLs are seen in our optical spectra, and a total 
carbon ORL abundance of 3.52 10 -3 by number with respect to 
hydrogen is derived. 

One magnesium line is seen in our spectrum - that of 
Mg II A4481. No ionisation correction scheme exists for magne- 
sium, but Mg 2+ occupies a very large ionisation potential interval, 
from ^15 eV to ~80eV, and so almost all the magnesium present 
will be in the form of Mg 2+ . The abundance derived for Mg 2+ is 
noteworthy. While all the second-row elements have ORL abun- 
dances which are much higher than those derived from CELs and 



exceed the solar abundance by factors of 10-30, the abundance de- 
rived for magnesium, using the model MC2's ICF (see Table 8 and 
Secions 4 and 5), is only a factor of 1.9 higher than the solar value. 
This is in line with the Mg abundances found by Barlow et al. 
(2003) for other nebulae showing large CNO ORL/CEL discrep- 
ancies, and strongly suggests that the large CNO ORL abundances 
are due to astrophysical effects such as H-deficient knots, and not to 
some unknown atomic process, since the latter might be expected 
to affect third-row as well as second-row elements. 

Argon, sulphur and chlorine abundances are derived from 
CELs only (see Table 6), with KB94 ionisation correction factors 
of 1.02, 2.73 and 2.83, respectively. 

The [Fe III] A4658 line is contaminated by the C IV A4658 
line (see Table 1). However, from the measured intensity of 
[Fe III] A4881 (Table 1), the second strongest line in the optical 
spectrum of [Fe III] , and the theoretical [Fe III] A4658/A488 1 inten- 
sity ratio of 2.0 at T c =12 kK and N c = 10 3 cm~ 3 (collision strength 
and A-value from Zhang 1996; Nahar 1996, respectively), we esti- 
mate that [Fe III] contributes 55% to the 4658 A blend, with C IV 
contributing the rest. This yields Fe 2+ /H+ = 3.82- 10~ 8 (see Ta- 
ble 6). Using an ICF of 155 from the model MC2, we derive a 
CEL-based total Fe abundance of 5.92 ■ 10~ 6 , by number with re- 
spect to hydrogen. 

We also calculated total abundances from the CELs using the 
ICFs from model MC2. Columns 3 and 5, in the upper half of Ta- 
ble 8, list the model ICFs and corresponding total abundances, ob- 
tained by applying the ICFs to the ionic abundances derived from 
the CEL analysis. Total elemental abundances obtained using the 
MC2 ICFs should be more accurate than those obtained by using 
the KB94 ICFs and are adopted in this work. The largest discrep- 
ancy between the KB94 ICFs and the model ones is found for nitro- 
gen. The CEL ionic abundances, together with the KB94 ICFs, give 
N/O = 0.20, which would make NGC 1501 a non-Type I PN ac- 
cording to both the criterion of Peimbert & Torres-Peimbert (1983, 
N/O > 0.5 for a Type I PN) and the criterion of KB94 (N/O > 0.85 
for a Milky Way Type I PN). However, the factor of four difference 
between the KB94 and MC2 ICFs for nitrogen turns NGC 1501 
into a Type I PN when using the MC2 ICFs, with (N/0) C el = 0.96. 

For the ORLs, when using MC2 ICFs we find N/O = 0.52, 
while for the KB94 ICFs we find N/O = 0.19. However, since CEL 
abundances are believed to provide a better representation of 'aver- 
age' nebular gas (see e.g. Pequignot et al. 2003), the classification 
of NGC 1501 as a Type I PN is preferred here. The elemental abun- 
dances derived from the CELs using MC2 ICFs are listed in Ta- 
ble 9, where they are compared with the ORL-based abundances, 
with average local disk PN abundances from KB94, and with so- 
lar abundances. For N, O and Ne, the ORL-based abundances are 
larger than the CEL-based abundances by factors of 7, 36 and 42, 
respectively. 



4 THREE-DIMENSIONAL PHOTOIONISATION 
MODELLING OF THE NEBULA 

4.1 The Computer Code and Model Setup 

The three-dimensional photoionisation code MOCASSIN- 
Version 1.10 (Ercolano et al. 003a) was used in order to 
construct a realistic model for NGC 1501. Sets of models were run, 
of which we selected two, MCI and MC2, which best reproduced 
the nebular emission spectrum. The initial models (MCI) were 
run on the HiPerSPACE facilities at University College London, 



The planetary nebula NGC 1501 1 3 



Table 9. Recommended elemental abundances in NGC 1501 from CELs 
and ORLs, in units such that log N(H) = 12.0 



Element 


ORLs 


CELs 


PN a 
Average 


Solar b 


He 


11.01 




11.06 


10.93 


C 


9.55 




8.74 


8.39 


N 


9.36 


8.52 


8.35 


7.92 


O 


10.09 


8.53 


8.68 


8.69 


Ne 


9.46 


7.84 


8.09 


8.08 


S 




7.52 


6.92 


7.33 


CI 




5.10 




5.50 


Ar 




6.11 


6.39 


6.40 


Fe 




6.77 




7.52 



a The overall average abundances of Galactic PNe (Kingsburgh & Barlow, 
1994); 

b Solar abundances from Grevesse & Sauval (1998, 1999), apart from C 
and O abundances which are from Allende Prieto et al. (2001, 2002). 



Table 10. Input parameters for the MOCASSIN photoionisation models. The 
abundances listed in this table are the final abundances used in models MCI 
and MC2 



Stellar and Nebular 
Parameters 



Nebular Abundances 
Model CEL 



T eB [kK] 


110 


He/H 


0.110 




logs (cgs) 


6.0 


C/H 


3.40(-4) 




U [L ] 


5000 


N/H 


3.161-4) 


3.28(-4) 


D [kpc] 


1.3 


O/H 


3.50(-4) 


3.42(-4) 


e 


0.423 


Ne/H 


6.35(-5) 


6.98(-5) 


^cubc [ ] 


38 


Mg/H 


3.80(-5) 








Si/H 


1.00(-5) 








S/H 


3.10(-5) 


3.31(-5) 






Cl/H 


1.20(-7) 


1.27 (-7) 






Ar/H 


1.10(-6) 


1.29(-6) 






Fe/H 


6.60(-6) 


5.92(-6) 



consisting of a SUN Microsystems V880 multiprocessor computer 
with 16 processors and 32Gb of memory. The later models (MC2) 
were run on a Beowulf cluster that belongs to the Institute of Com- 
putational Mathematics and Scientific/Engineering Computing of 
the Chinese Academy of Sciences, consisting of 120 550-MHz 
Pentium III processors with 1Gb of memory each. A grid of 71 3 
spatial points was used for our final models, corresponding to 
357911 cubic cells of length 1405 AU each (assuming a distance 
of 1300 pc). The only difference between the MCI and MC2 
models is that the latter uses test values (see Section 5.1.1) for the 
unknown low-temperature dielectronic recombination coefficients 
of third-row elements (S, CI and Ar), while the former sets all 
unavailable rates to zero. The nebular and stellar input parameters, 
as well as the final model abundances for MCI and MC2, are listed 
in Table 10. 



4.2 The Density Distribution 

The three-dimensional density distribution grid used for the mod- 
elling was determined for by Ragazzoni et al. (2001) from long-slit 
echellograms, reduced according to the methodology developed by 



Sabbadin et al. (000a,b). As the density distribution is, in general, 
a difficult input parameter to constrain, given the plurality of 
acceptable solutions, the existence, in this case, of a model density 
directly derived from the observations puts us in a favourable 
position, allowing a potentially higher degree of accuracy in the 
determination of the other nebular and central star parameters. The 
data cube made available to us by F. Sabbadin consisted of a 200 3 
grid of N e ; these were mapped to the 71 3 MOCASSIN grid and 
converted to gas densities by iteratively using the local ionisation 
structure calculated at each grid cell to calculate the conversion 
factors. 

Figures 4 and 5 of Ragazzoni et al. (2001) show iso-density 
plots of the 7V e distribution they derived at two different threshold 
values (N e = 300 cm -3 and N c = 900 cm -3 ), as seen from 12 dif- 
ferent directions, while the optical appearance of the rebuilt nebula 
as seen from the same 12 directions is shown in their Figure 6. It is 
worth noting at this point that the original N c cube contained val- 
ues of TVe(Ha), where N c (Ha) = N e (S II) ■ e 5 , e being the local 
filling factor. Unfortunately, the value of e could not be determined 
from Ragazzoni et al.'s observations. The value of 0.423 used in 
the MOCASSIN models (and reported in Table 10) was chosen in or- 
der to reproduce /(H/3) = 5.26 10~ n ergs -1 cm -2 , the integrated 
H(/3) flux from Cahn et al. (1992), dereddened using c(H/3)= 1.0 
(see Section 3.3). This was made possible by the fact that this large, 
old, low-electron density PN is optically thin in H I, implying that 
the H/3 emission of the nebula is not dependent on the luminos- 
ity of the central star, but represents a measure of the ionised mass. 

From this it is clear that I(Ha) oc -gr, where D is the distance from 
the nebula. As is true for many astronomical objects, the distance 
to NGC 1501 is quite poorly constrained. At least 15 estimates are 
available in the literature, including statistical and individual dis- 
tance estimates, with values ranging from 0.9 kpc (Amnuel et al. 
1984) to 2.0 kpc (Acker 1978). The value of 1.30 kpc adopted in 
this work, was chosen for reasons of consistency with the value 
adopted by Ragazzoni et al. (2001) for their calculation of the three- 
dimensional ionisation structure from observations. We finally note 
that, although S + may not be a faithful representation of the ionised 
gas in this high-excitation nebula, Ragazzoni et al.'s result is cer- 
tainly an improvement from the assumption of a spherically sym- 
metric analytical density law. 



4.3 Elemental Abundances 

A homogeneous elemental abundance distribution was used for our 
models given that, although the ORL-CEL ADFs are high in this 
nebula, which may indicate the presence of cold ionised hydrogen- 
deficient material (e.g. Liu et al. 2000), we do not have any direct 
observational evidence for the presence of structures containing 
plasma with different elemental abundances. It is, however, pos- 
sible that such structures might exist inside the nebular shell, but 
with sizes falling below the spatial resolution of the instruments 
currently available. With the HST, a spatial resolution of 0.045"can 
be achieved in the optical - this corresponds to 53 AU at 1300 pc. 
HST/WFPC2 images of NGC 1501 have been taken as part of GO 
program 6119. However, the observations were aimed at detecting 
companions to central stars, and the short exposures necessary to 
obtain good stellar images yield very poor signal-to-noise ratios 
for the nebula itself. Also, the observations were taken in broad- 
band filters F555W and F814W, which contain many emission lines 
and are thus impossible to use for quantitative nebular analysis. 
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As will be discussed in more detail in Section 6, the presence of 
small, hydrogen-deficient knots could provide an explanation for 
the aforementioned ORL-CEL discrepancies as well as other spec- 
troscopic peculiarities. 

The abundances derived empirically from CELs (see Sec- 
tion 3.5) were used as starting values for our modelling; these were 
then modified in order to obtain a better fit to the CEL spectrum. 



4.4 Central Star Parameters 

The effective temperature and gravity of the central star were var- 
ied in order to provide an ionising spectrum that could reproduce 
the ionisation structure implied by the observed nebular spectrum, 
under the constraints of our adopted density distribution, filling fac- 
tor and distance. The final values adopted are T e g = 110kK and 
log g = 6.0 (cgs). This is well within the range of effective temper- 
ature determinations available in the literature, which range from 
Tzan S tra(He II) = 81kK (Stanghellini et al. 1994) to the value of 
135kK obtained by Koesterke & Hamann (1997) by means of stel- 
lar atmosphere modelling. It is worth noting at this point that, 
within a certain range of gravities, the value of log g has a negli- 
gible influence on the ionising flux level of these models. 

The hydrogen-deficient non-LTE stellar atmospheres used in 
this work were calculated using the Tubingen NLTE Model At- 
mosphere Package (Rauch 2003), for He:C::N:0 = 33:50:2:15 by 
mass, i.e. abundances close to those derived for the central star 
of NGC 1501 and for PG 1159 stars (see Section 3.4). These 
models are available for download from Dr. T. Rauch's website 
(http : // astro . uni-tuebingen . de/~rauch). The best 
results for the nebular photoionisation modelling were obtained 
with a central star effective temperature of 110 kK (as mentioned 
above) and a stellar luminosity of L, = 5000 L Q , which is close to 
the 5500 Lq estimate of Koesterke & Hamann (1997). This model 
was found to be consistent with the level of the stellar UV contin- 
uum between 1300-1800 A in the dereddened IUE spectrum (see 
Section 2.2). 



4.5 The Ionising Spectrum 

As discussed in Section 3.2, evidence for a strong stellar wind is 
observed in the spectrum of the central star of NGC 1501. This 
poses the question of whether a model atmosphere calculated with 
the plane parallel Tubingen code is a justifiable choice. We com- 
pared it with a model atmosphere calculated with the spherically 
symmetric code CMFgen (Hillier & Miller 1998), for wind pa- 
rameters close to those of the central star of NGC 1501 and abun- 
dances He:C:N:0 = 54:36:1:8 by mass, similar to those derived by 
Koesteke & Hamman (1997; He:C:0 = 50:35: 15), with a Tubingen 
model for the PG1159-type abundances given above; both models 
had T c ff = 135kK. Despite the different abundances, the variations 
in the flux levels between the two models were found to be less 
than a factor of two, for wavelengths shortward of the Lyman limit, 
which is within the error bars imposed by the uncertainties of the 
stellar abundance determination. We also compared our Tubingen 
model atmosphere at HOkK with one calculated by L. Koesteke 
(private communication) using a spherically symmetric wind model 
for the same temperature and PG1 159-type abundances, and found 
that the flux levels predicted by the two models agreed to within 
20%. 



5 NEBULAR MODELLING: STRATEGY AND RESULTS 

The results of our photoionisation modelling and a comparison 
with observed values are presented in this Section. Our preliminary 
modelling was aimed at constraining the basic stellar and nebular 
parameters. Bearing in mind that the empirical analysis of the ORL 
spectrum of NGC 1501 yields larger abundances than the CEL 
analysis, and that this could be interpreted as indicating the exis- 
tence of a separate component in the nebular gas, our efforts aimed 
to match the CEL spectrum only, as this should be representative of 
the dominant mass component of the nebular gas (Liu et al. 2000; 
Pequignot et al. 2003). 



5.1 The first model (MCI) 

The initial model will be referred to in this paper as MCI. The 
CEL intensities predicted by MCI are compared to the observed 
values in Table 11. Table 12 lists the predicted and observed nebu- 
lar diagnostic line ratios. The results obtained from this first round 
of modelling, of which only the best-fitting model is shown in 
the table, were quite surprising. A comparison of the CEL inten- 
sities predicted by the model with the observations indicates that 
the model is too highly ionised. However this is at odds with what 
is implied by the predicted He I and He II recombination lines, 
which show that the degree of ionisation of He is actually slightly 
lower than indicated by observations (see Table 13). The problem 
is evident from examination of the values listed in Table 14, con- 
taining the X l+1 /X +l model ionic fractions (columns 2 and 3) and 
those derived from the empirical analysis (column 4). Changes of 
the central star ionising spectrum, aimed to adjust the ionisation 
degree of O, S, CI and Ar in the nebula, would necessarily also 
affect the ionisation structure of He, thus causing a larger discrep- 
ancy with the observations. This was confirmed by a comparison of 
models run with several T ff-log g combinations within the range 
T eS = 100kK-140kK and log g = 5.0-8.0 [cgs]. 



5.7.7 The problem of data incompleteness for the 

low-temperature dielectronic recombination process 

The over-ionisation problem of S, CI, and Ar could be interpreted 
as evidence for the effects of low-temperature dielectronic recom- 
bination, a process which can in magnitude be comparable (or often 
larger) than its radiative counterpart, but for which no reliable rates 
have yet been calculated for third-row elements. The uncertainty 
of the atomic data available for some astrophysically abundant ele- 
ments is faced by all photoionisation codes. Ali et al. (1991) argued 
that, for a given ion X + \ a better estimate than zero is obtained by 
taking the average of the rate coefficients for C + \ +I and N + \ 
given that the coefficients seem to follow a certain trend with ioni- 
sation stage. Dudziak et al. (2000), use coefficients empirically cal- 
ibrated from a new unpublished model of the PN NGC 7027. In 
this work, we do not attempt to calculate empirical coefficients, as 
NGC 1501, with its complicated nebular spectrum and WR central 
star, cannot be considered a typical PN. Self-consistent calculations 
of dielectronic and radiative rates of some third-row elements are 
currently being carried out (Storey, Ercolano and Badnell, in prepa- 
ration), but in this paper we will limit ourselves to seeking empir- 
ical upper limits to the unknown dielectronic rates and use these 
coefficients as test values for our nebular modelling. 

For a pocket of gas in ionisation equilibrium, the abundance 
ratio of two successive ionic stages is given by 
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Table 12. Calculated and observed plasma diagnostic ratios 



Ion 


Lines 


IP(eV) 


Observed 


MCI 


MC2 


[S II] 


7(A6731)//(A6716) 


10.36 


1.10 


1.12 


1.12 


[0 ii] 


i"(A3729)/i"(A3726) 


13.62 


0.83 


0.86 


0.86 


[On] 


7"(A3727)//(A7325) 


13.62 


16.85 


26.84 


26.11 


[CI in] 


7(A5537)//(A5517) 


23.81 


0.95 


0.84 


0.84 


[Ar IV] 


7(A4740)/7(A4711) 


40.74 


0.75 


0.78 


0.77 


[Nil] 


/(A6548 + A6584)//(A5754) 


14.53 


52.38 


68.67 


68.60 


[0 III] 


/(A4959 + A5007)//(A4363) 


35.12 


148.7 


138.7 


138.4 


He I 


7(A5876)//(A4471) 


24.59 


3.03 


2.72 


2.72 


He I 


/(A6678)//(A4471) 


24.59 


0.86 


0.76 


0.76 



Table 11. Collisionally excited emission line intensities predicted by 
MOCASSIN, compared to observed values. The line intensities are given rel- 
ative to 7 (H/3) on a scale where i"(H/3) = 100. The intensity of H/3 is given 
in absolute units [erg cm~ 2 s -1 ] and is compared to the value given by 
(Cahn et al. 1992), dereddened by c(H/3) = 1.0. 



Table 13. Helium recombination line intensities predicted by MOCASSIN, 
compared to observed values. The line intensities are given relative to 
I (H/3) on a scale where 7(H/3) = 100. The intensity of H/3 is given in ab- 
solute units [erg cm~ 2 s~ x ] and is compared to the value given by (Cahn 
et al. 1992), dereddened by c(H/3) = 1.0. 



X 



i+l 



Line 




AA CI 


Observed 


Log[I(H/?)] 


-10.27 


-10.27 


-10.28 


H/3 A4861 


100 


100 


100 


[N II] A5755 


0.20 


0.20 


0.25 


[N II] A6548 


3.44 


3.45 


3.16 


[N II] A6584 


10.14 


10.01 


9.94 


[O II] A3727 


5.44 


10.7 


12.16 


[O II] A3729 


4.69 


9.23 


10.08 


[O II] AA7318.9 


0.21 


0.41 


0.73 


[O II] AA7330.0 


0.17 


0.33 


0.59 


[O III] A4363 


12.81 


12.78 


10.38 


[O in] A4932 


0.19 


0.19 


0.35 


[O III] A4959 


458 


456 


394 


[O III] A5007 


1319 


1312 


1151 


[Ne III] A3869 


102 


103 


101 


[Ne III] A3967 


31.6 


31.7 


28.10 


[S II] A4069 


0.069 


0.28 


1.49 


[S II] A4076 


0.023 


0.10 


0.32 


[S II] A6717 


0.43 


1.76 


1.77 


[S ii] A6731 


0.48 


1.98 


1.95 


[S III] A6312* 


1.48 


1.48 


1.48 


[CI III] A5518 


0.11 


0.84 


0.76 


[CI III] A5538 


0.090 


0.71 


0.72 


[CI IV] A7531 


0.56 


0.26 


0.26 


[CI IV] A8046 


1.31 


0.61 


0.54 


[Ar III] A7136 


1.30 


13.08 


14.34 


[Ar III] A7752 


0.31 


3.16 


3.52 


[Ariv] A4712 


9.02 


2.74 


2.66 


[Ar IV] A4741 


7.01 


2.12 


2.00 



* Blend with Hell A6310 



X i 



N x 

Jv phot 

N e ■ (a R (X+ l ,T c ) + a D (X+ l ,T c )) 



(2) 



Line 


MCI 


MC2 


Observed 


Log[I(H/3)] 


-10.27 


-10.27 


-10.28 


H/3 A4861 


100 


100 


100 


He I A4471 


4.20 


4.20 


3.76 


He I A4922 


1.11 


1.11 


1.07 


He I A5876 


11.42 


11.43 


11.40 


He I A6678 


3.21 


3.21 


3.12 


He I A7065 


2.16 


2.16 


2.51 


He II A4686 


34.7 


34.6 


40.30 



Table 14. Predicted and observed ionic ratios 



Xi+l/X* 


MCI 


MC2 


obs 


MCl/obs 


MC2/obs 


2 +/0+ 


116 


58.9 


46.5 


2.49 


1.27 


s 2 +/s+ 


46.8 


27.8 


28.3 


1.65 


0.98 


C1 3 +/C1 2 + 


10.4 


0.61 


0.60 


17.4 


1.02 


Ar 3 +/Ar 2 + 


12.0 


0.35 


0.33 


36.2 


1.06 




MCI 


MC2 


obs 


MCl/obs 


MC2/obs 


He 2 +/He+ 


0.36 


0.36 


0.41 


0.88 


0.88 



Table 15. Estimated upper limits to low-temperature dielectronic rates 







a D (X\10 4 K) 






[cnv^s- 1 ] 


s+ 


3.14 


5.66(-12) 


Cl 2 + 


17.5 


5.81(-11) 


Ar 2 + 


42.5 


13.7(-11) 



a f(Xi) = a D (X\ W 4 K)/a R (X\ 10 4 K) 



where N pbot is the total number of photoionisation events and 
ol r {X +1 ,T e ) and a D (X +l ,T e ) are the temperature-dependent 
rate coefficients for, respectively, radiative and dielectronic recom- 



16 Ercolano et al. 



Table 16. Low-temperature dielectronic rates for CNO ions (Nussbaumer 
& Storey 1984) 



X 1 


f(x i )a 


[crn^- 1 ] 


c+ 


2.4 


6.06(-12) 


C 2+ 


2.7 


1.36(-11) 


N+ 


.86 


2.03(-12) 


N 2 + 


4.1 


2.16(-11) 


0+ 


.83 


1.66(-12) 


2 + 


2.0 


1.14(-11) 



"ftX*) = a D (X\ W 4 K)/a R (X\ 10 4 K) 



bination to X + \ Note that we are assuming that the contribution 
due to charge exchange processes is negligible in this case, given 
that the nebula is optically thin in H I. This implies that we can use 
the ratio of model to observed X'+VX 1 values to obtain an estimate 

of a D {X+\T c ): 

(( y"* + 1 / Y'-\-i\m.od \ 
\ x J /x+ : )ob s -lj -<*n{X + \T*) 0) 

Unfortunately, our observations did not yield all the X* +1 /X +I 
ratios necessary to derive a consistent set of estimates for all 
the relevant ionic stages of S, CI and Ar, but only for S + , 
Ar 2+ and Cl 2+ . We define a quantity f(X +i ) such that the 
dielectronic rates for each ionic stage can be expressed as 
a D (X+\10 4 K) = /(X +i )a J? .(X+\10 4 K), where initial guesses 
for f(X +% ) were taken from Equation 3 and then adjusted itera- 
tively with the modelling. From the above, f(X +l ) is defined as 
the ratio of the low-temperature dielectronic recombination coef- 
ficient to the radiative recombination coefficient of the ion X 1 at 
T c = 10 4 K. The use of a temperature independent coefficient such 
as f(X +t ) is equivalent to saying that the dielectronic rates follow 
the same temperature dependence as the radiative rates for the same 
ion, which is, of course, not true. Nevertheless, this is the best we 
can do in this case, as we are not attempting to calculate precise 
rates, a task that can only be fulfilled by laboratory measurements 
or detailed atomic modelling of each individual ion. Furthermore, 
the nebular gas T 's do not show large deviations from a mean of 
1 1 100 K, and therefore the error introduced by the temperature de- 
pendence will be much smaller than the margin of error implied by 
our empirical approach. The final f(X +l ) values used are given in 
Table 15, together with the corresponding values of ac(I + ', T c ) 
for T e = 10 4 K. These values are only to be considered as upper lim- 
its, since other effects may be at play in complicating the ionisation 
structure of this object (see Section 6). As a comparison, the low- 
temperature dielectronic recombination coefficients and the corre- 
sponding / values are listed in Table 16 for C, N and O ions (Nuss- 
baumer & Storey 1984) at lOkK. Using the Ali et al. (1991) method 
method described above with the values listed in Table 16, one ob- 
tains estimates of 3.25 ■ 10~ 12 cm 3 s _1 and 1.55 • 10~ n cm 3 s _1 for 
the singly- and doubly-ionised cases at T c = lOkK, lower than the 
values listed in Table 15, necessary to resolve the discrepancy with 
the observations. 



5.2 Models including estimated low-temperature dielectronic 
recombination rates for third-row elements (MC2) 

5.2.1 The CEL spectrum 

The CEL intensities predicted by the best-fitting model from a sec- 
ond set of models (MC2), which use the estimated upper limits to 
the dielectronic recombination rates given in Table 15, are given 
in column 3 of Table 1 1 . This model provides a much better fit to 
the overall observed emission line spectrum than model MCI, but 
some problems are still present, as discussed in Section 6. 

The ionisation degree of O also appeared to be overestimated 
by model MCI (see Table 14). The dielectronic rates included 
in the atomic database of MOCASSIN, were obtained by Nuss- 
baumer & Storey (1984) who performed the calculations within 
the LS coupling scheme. As discussed in Section 4.3 of the Nuss- 
baumer & Storey (1984) paper, using intermediate coupling (IC) 
(e.g. Beigman & Chichkov 1980) could result in higher rates as 
more states are taken into account (although not all of these would 
be permitted to autoionise). With this in mind, we increased the 
dielectronic rates in MC2 by a factor of three, which returned a 
better fit of the oxygen CEL spectrum. However, as discussed in 
Section 6, the large corrections to the atomic data required by MC2 
are questionable and highlight the limitations of this chemically 
homogeneous model. 

The [O II] A7319 and A7330 doublets and the [S II] A4069 and 
A4076 lines are still underestimated by model MC2. Recombina- 
tion processes can contribute to the intensity of these lines, and, 
whilst effective recombination coefficients to the metastable lev- 
els of [O II] are known (Liu et al. 2000), the same is not true for 
[S II]. We have included the recombination term to the statistical 
equilibrium matrix in our photionisation model, in order to self- 
consistently calculate the populations of the metastable levels of 
[O II], by taking into account recombination and collisional popu- 
lation and depopulation. However, at the temperatures and densities 
of the nebular gas, we found recombination contributions to be neg- 
ligible. We therefore conclude that the discrepancy between model 
and observed intensities of these line must be attributed to recom- 
bination processes occurring in a separate, colder and/or denser gas 
phase (see Section 6). 



5.2.2 The computed ionic and thermal structure of the nebula 

The nebular-averaged fractional ionic abundances calculated from 
the MC2 model are listed in Table 17. Hydrogen and helium are 
both fully ionised and we see that significant fractions of the heavy 
elements are in higher ionisation stages. We notice that the N/N + 
ratio is nearly a factor of three larger than the 0/0 + ratio, contrary 
to what is generally assumed in ionisation correction schemes, in- 
cluding the one used in this work, where the two ionic fractions are 
taken to be the same. This has already been noticed by Rubin et al. 
(1988) and should be taken into account when considering uncer- 
tainties in empirical abundance analyses for optically thin PNe such 
as NGC 1501 (see also Section 3.7). 

We did not find large temperature fluctuations in the computed 
3D temperature distribution, obtaining t 2 (0 2+ ) = 0.0021 from the 
model. This is not surprising, given that large density fluctuations 
are not present in our density distribution grid. 
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Table 17. Nebular-averaged fractional ionic abundances, obtained from model MC2. 



Ion 



Element 


I 


II 


III 


IV 


V 


VI 


VII 


H 


1.27(-3) 


9.99(-l) 












He 


7.43(-4) 


7.34(-l) 


2.65(-l) 










C 


4.60(-6) 


8.89(-3) 


4.09(-l) 


5.67(-l) 


1.49(-2) 






N 


9.95(-7) 


3.97(-3) 


3.91(-1) 


6.02(-l) 


3.11 (-3) 






O 


1.24(-5) 


1.47(-2) 


8.66(-l) 


1.18(-1) 


1.31 (-3) 






Ne 


9.37(-7) 


5.03(-3) 


9.42(-l) 


5.30(-2) 


3.36(-5) 






Mg 


2.98(-3) 


3.11(-2) 


9.15(-1) 


5.12(-2) 


1.93(-6) 






Si 


5.40(-5) 


3.77(-2) 


2.64(-l) 


4.18(-1) 


2.80(-l) 






S 


3.36(-7) 


2.42(-3) 


6.72(-2) 


4.87(-l) 


4.39(-l) 


4.41(-3) 


5.89(-6) 


CI 


1.13(-6) 


4.44(-3) 


5.51(-1) 


3.57(-l) 


1.05(-1) 


2.94(-3) 


1.12(-7) 


Ar 


1.18(-7) 


9.40(-4) 


7.14(-1) 


2.53(-l) 


2.90(-2) 


2.63(-3) 


9.37(-6) 


Fe 


4.05(-5) 


6.72(-4) 


6.44(-3) 


6.47(-l) 


3.13(-1) 


3.22(-2) 


2.15(-5) 



6 THE USUAL SUSPECTS: t 2 VERSUS 

ENHANCED-DENSITY, METAL-RICH KNOTS 

The introduction of the test values for the low-temperature dielec- 
tronic recombination coefficients of S, Ar and CI certainly helped 
to achieve a better fit to the observed CEL spectrum of NGC 1501. 
However a number of problems still remain unsolved. In particular: 

1. The ORL-CEL discrepancy. MC2 fails to reproduce the ob- 
served ORL spectrum. The abundances derived empirically from 
the ORL spectrum are consistently larger than those derived from 
the CEL analysis, while the T c derived from the observed He I re- 
combination lines is cooler than the temperature estimates obtained 
from CELs (Table 5). 

2. The large dielectronic rate value for argon. With regards to the 
Ar 3+ /Ar 2+ ionic ratios, large low-temperature dielectronic recom- 
bination rates are required in order to resolve the discrepancy be- 
tween the observations and the photoionisation models. Although 
strongly dependent on the atomic structure of each individual ion, 
atomic data currently available for a number of doubly ionised ions 
show low-temperature dielectronic rates to be at most ten times 
larger than the corresponding radiative rates (with average values 
being in the region of 3 or 4 times the radiative rates). On the other 
hand, recent work by Morisset et al. (2004) also pointed out the 
need for large dielectronic rates for Ar + in order to reconcile ISO 
infrared observations and models for forty H II regions in their sam- 
ple. 

Temperature fluctuations have been proposed within the t 2 
framework of Peimbert (1967, 1971, and later references) by 
Ruiz et al. (e.g. 2003) and Peimbert et al. (2004) as a solution to 
the ORL-CEL abundance discrepancy. In the case of NGC 1501 
our three-dimensional models showed no evidence for t 2 fluctu- 
ations larger than 0.012. We also calculated empirical values for 
T ([O III]) and t 2 ([0 III]), following the formulations given by Pe- 
imbert (1967, 1971) and, more recently Peimbert et al. (2004). First 
of all we calculated that a value of T)([0 III]) = 3750 K is required 
to obtain the same abundances from the [O III] CELs and the O II 
ORLs, implying 43% temperature fluctuations (t 2 ([0 III]) = 0.184). 
However, when we derived values for the same quantities using the 
Balmer jump and the [O III] temperatures from Table 5, we found a 
mean temperature of T ([O III]) = 10010 K and t 2 ([0 III]) = 0.036, 
which correspond to temperature fluctuations of only 19%. We also 



calculated To and t 2 values using the T c derived from the He I 
recombination lines (see Table 5). The He I A5876/A4471 T c of 
4100 K implies that T = 5350 K and t 2 = 0.154, while the mean 
He I T c of 5 125 K implies that T = 6300 K and t 2 = 0.133. These 
values of To and t 2 would increase the 2+ abundance measured 
from CELs by factors of 10.6, 3.3 and 6.5, respectively, still a long 
way short of resolving the factor of 30 discrepancy. The observa- 
tional results and the computed 3D temperature structure indicate 
that something other than temperature fluctuations is causing the 
large ORL-CEL ADFs we found for NGC 1501. 

The existence of metal-rich knots could provide a solution to 
the problems mentioned above. While such knots have been in- 
voked recently to help with the solution of the ORL-CEL discrep- 
ancy problem in PNe (e.g. Liu, 2000), and bi-abundance photoioni- 
sation models have been constructed for some PNe (Pequignot et al. 
2003), it is further proposed here that the lower ionisation degree 
of heavy elements compared to helium shown by NGC 1501 can 
also be explained by the presence of small metal-rich knots mixed 
with the nebular gas. The knots could be formed by processed ma- 
terial transported out in the outflow from the AGB precursor of the 
H-deficient WR central star. 

Low-ionisation species would be abundant in the knots and in 
the shadow regions (tails) behind them, thanks to softening of the 
radiation field due to the screening effect of material that is opti- 
cally thick to stellar UV radiation. However, due to the high metal- 
licity and cooling rates in the material, the temperatures inside the 
knots are expected to be much too low for collisionally excited line 
radiation to be produced to a significant degree, while essentially all 
the ORL emission is expected to be produced in these regions (see 
e.g. the 3-D photoionisation model of the H-poor knots in Abell 30 
by Ercolano et al. 003b). 

The situation in the shadowed tails of knots is completely dif- 
ferent. Since these are composed of material with normal nebular 
abundances, they are not as efficiently cooled as the knots them- 
selves, implying that the emission of collisionally excited lines 
from lower ionisation species would be expected. In this regard it 
is clear that the problem of over-ionisation of oxygen in the CEL- 
emitting region for models MCI and MC2 could be eased by the 
introduction of such knots, which could also provide some of the 
missing emission from S + , Ar + and CI + , with the implication 
that lower dielectronic rates than those listed in Table 15 would be 
sufficient to match the CEL spectrum. 
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7 CONCLUSIONS 

We have presented in this paper new, deep optical spectra of the 
PN NGC 1501. The empirical analysis led to the discovery of high 
ORL-CEL ADFs for 2+ and Ne 2+ , with values of 32 and 33, re- 
spectively. The agreement between the ADFs for O and Ne lends 
further weight to the suggestion by Liu et al. (2000) that, while 
large discrepancies may exist between ORL and CEL abundances 
relative to H, elemental ratios such as Ne/O are often the same 
whether derived from CELs or ORLs. 

Unfortunately no temperature for the ORL-emitting region 
could be determined from our observations, since the O II A4089 
line (see Wesson et al. 2003; Tsamis et al. 2004) was not detected 
in our observations. However, a large temperature discrepancy ex- 
ists between the T c implied by the He I recombination line ratios 
(on average ~5.1 kK) and the [O III] temperature (11.4 kK). Given 
that we expect T c (ORL) < T c (He I) < T c (CEL) (Liu, 2003), 
the ORL emitting regions probably have temperatures lower than 
4.5 kK, consistent with the results of Wesson et al. (2003) for knot 
J3 of Abell 30, where T (ORL) = 2.5 kK. T e (He I) > T c (ORL) 
is a logical consequence of He I recombination lines being pro- 
duced in both the ORL and the CEL-emitting regions and, there- 
fore, the T e 's derived from them are not representative of either 
region. 

Our initial models showed metals to be too highly ionised 
compared to observations, except for helium. We suspect that 
two effects are at play in complicating the emission line spec- 
trum of NGC 1501. First of all, we note the incompleteness of 
low-temperature dielectronic recombination coefficients for third- 
row elements and the possible underestimation of qd(0 + , T). 
We have obtained upper limits for a D (S+, 10 4 K), a c (Cl 2+ , 
10 4 K) and a c (Ar 2+ , 10 4 K) of 5.66-10" 12 , 5.8110" 11 and 
13.7-10 - 1 cm s~ , respectively. These correspond to low- 
temperature dielectronic recombination coefficients that are factors 
of 3.1, 17.5 and 42.5 higher than the corresponding radiative re- 
combination rates for S + , Cl 2+ and Ar 2+ , respectively. We also 
found that increasing q?d(0 + , T) by a factor of three yields a better 
fit to the [O II] spectrum. The very large rates required for Ar and 
CI make it improbable that the discrepancies between the models 
and observations can be resolved solely in terms of these coeffi- 
cients. This, together with the large ORL-CEL ADFs derived for 
NGC 1501, points to the existence of a second component, consist- 
ing of cold, metal-rich ORL-emitting knots. It seems that the solu- 
tion to the over-ionisation problem may lie partly with the softening 
of the ionising radiation field, due to absorption of UV photons by 
such knots and partly with the lack/uncertainty of low-temperature 
dielectronic rates. Clearly the magnitudes of the dielectronic rates 
required to fit the spectrum could be reduced by an amount de- 
pending on the fraction of the emitting gas comprised by the knots, 
as well as their size and density, which will affect the size of the 
shadowed regions. 

We also note that NGC 1501 is yet another PN with a 
hydrogen-deficient central star which shows high ADFs. To our 
knowledge, NGC 1501, with its ADF of ~32, has the third highest 
ADF of any nebula, after Abell 30 (ADF~700; Wesson et al. 
2003) and Hf 2-2 (ADF-84; Liu 2003). Hydrogen-deficient knots 
are directly observed in the former (lacoby et al., 1979; Hazard 
et al., 1980) and three-dimensional photoionisation modelling 
(Ercolano et al., 2003b) has shown that in its case cold, metal-rich 
material can account for the ORL emission observed. In the case 
of Abell 30, a popular view with regards to the origins of such 
knots is the born-again scenario (Iben & Renzini 1983), whereby 



hydrogen-deficient material would have been stripped off the 
stellar surface during a late thermal pulse event. 

There are, however, inconsistencies between the predictions 
of the born-again scenario and observations of some of the nebulae 
it applies to (including Abell 30). Most importantly, evolutionary 
models predict a late thermal pulse event to produce a carbon-rich 
stellar surface abundance (Herwig et al. 1999), whereas in the 
case of Abell 30 Wesson et al. (2003) measured (C/O) or i = 0.8; 
similarly, for NGC 1501 we measured (C/O) or i = 0.31. A rec- 
onciliation of theory and observations could be obtained if the 
ORLs were produced in a gas component consisting of carbon-rich 
material expelled from the stellar surface, which has already 
mixed somewhat with surrounding oxygen-rich nebular gas. In 
the case of a dense knot, one would expect the ionising radiation 
to be able to penetrate only a thin outer layer, the same boundary 
region would also be the most likely to be subject to mixing. If 
the above is true, one would then expect the measured nebular 
(C/O) ri to be larger than (C/0) cc i, which is indeed confirmed 
by the data available in the literature. In particular we looked 
at a ORL-CEL study of 16 PNe by Tsamis et al. (2004), for 
which ten PNe had (C/0) O ri/(C/0) co i >1, two lacked a value for 
(C/O) or i, two lacked a value for (C/0) cc i, and only for two PNe 
was (C/0)ori/(C/0) cc i <1. A similar trend is also confirmed by the 
work of Wesson (2004). 

If the origin of the hydrogen-deficient knots can really be traced to 
material ejected from the surface of the AGB precursors of the cen- 
tral stars of born-again nebulae, one might also expect a correlation 
between hydrogen-deficient central stars and nebulae showing 
large ORL-CEL abundance discrepancies. Such a correlation did 
not emerge from investigation of the measured ORL-CEL ADFs 
(Liu et al. 1995, 2000, 001b; Liu 2003; De Marco & Barlow 2001; 
Tsamis et al. 2004, Wesson 2004) and central star classification 
(e.g. De Marco & Barlow 2001) available in the literature. On 
the one hand, the central star spectral type is not always known 
and we cannot therefore exclude that some high- ADF central stars 
are not actually H-deficient. For instance, the spectral type of the 
central star of the PN Hf 2-2 (ADF=84; Liu 2003) is unknown. On 
the other hand, the IUE spectra of the central star of PN DdDm 1 
(ADF=11.8 Wesson 2004) has no emission lines of C IV A1550 
or He II A 1640, which are expected in WR spectra. The jury is 
therefore still out as to whether the concept of H-deficient knots 
as an explanation for large ADFs can be reconciled with current 
post-AGB stellar evolutionary theories. We can however be certain 
that the high ADFs in nebulae such as NGC 1501 cannot be 
explained by nebular temperature fluctuations of the classically 
conceived type. 
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